Preprint typeset using WTgfi. style cmulatcapj v. 04/21/05 



STOCHASTIC NATURE OF GRAVITATIONAL WAVES FROM SUPERNOVA EXPLOSIONS WITH 

STANDING ACCRETION SHOCK INSTABILITY 

Kei Kotake 1 , Wakana Iwakami 2 , Naofumi Ohnishi 2 , and Shoichi Yamada 3,4 



ABSTRACT 

We study properties of gravitational waves based on the three-dimensional simulations, which 
demonstrate the neutrino-driven explosions aided by the standing accretion shock instability (SASI). 
Pushed by evidence supporting slow rotation prior to core-collapse, we focus on the asphericities in neu- 
trino emissions and matter motions outside the protoneutron star. By performing a ray-tracing calcu- 
lation in 3D, we estimate accurately the gravitational waveforms from anisotropic neutrino emissions. 
In contrast to the previous work assuming axisymmetry, we find that the gravitational waveforms 
vary much more stochastically because the explosion anisotropies depend sensitively on the growth of 
the SASI which develops chaotically in all directions. Our results show that the gravitational-wave 
spectrum has its peak near ~ 100 Hz, reflecting the SASI-induced matter overturns of ~ O(10) ms. 
We point out that the detection of such signals, possibly visible to the LIGO-class detectors for a 
Galactic supernova, could be an important probe into the long- veiled explosion mechanism. 
Subject headings: supernovae: general — gravitational waves — neutrinos — hydrodynamics — in- 
stability 



1. INTRODUCTION 

The gravitational- wave astronomy is now becoming re- 
ality. In fact, significant progress has been made on 
the gravitational wave detectors, such as LIGO, VIGRO, 
GEO600, TAMA300, and AIGO with t heir international 
netwo rk of the observatories (e.g., iHough fe Rowanl 
(|2005f ) for a review). For the detectors, core-collapse su- 
pernovae especially in our Galaxy, have been supposed as 
one of the most pl ausible source s of gr a vitati o nal w aves 
(see, for example, IKotake et all (|2006f l: lOttl (|2009t ) for 
recent reviews). 

Traditionally, most of the model calculations of grav- 
itational waves (GWs) have focused on the bounce sig- 
nals in th e context of r ota tional (e. g.. IZwerger fe Mulled 
19971) ; IKotake et ail (l2003f): IShibata fc Sekiguchil 
2001; lOtt et alJ (120041 ): iDimmelmeier et all (|2007f ); 
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and refe rences therein) and 
core-coll apse (IKotake et ail 12004 
lObergaulinger et alJ 120061 ; ICerda Duran et all I2007T ). 
However recent stellar evolution calculations suggest 
that such rapid rotation assumed in most of the 
previous studies, albeit attractin g much attention 
currently as a relevance to collapsar fH irschi et al J 120051 : 
IWooslev fc Heg er 2006), is not canonical for progenitors 
of core-collapse supernovae with ne utron star formations 
dHeger et al.ll2005TlOtt et al.ll2006T ) 

In the case of the slowly rotating supernova cores, two 
other ingredients are expected to be important in the 
much later phases after core bounce, namely convective 
motions and anisotropic neutrino emissions. Thus far, 
various physical ingredients for producing asphericities 
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and the resulting GWs in the postbounce phase, have 
been studied s uch as the roles of pre-co l lapse density in- 
homogeneities (jBurrows fc Ha ves 1996; M filler fc Jarikal 
119971: iFrverl l2004h. moderate rotation of the iron core 
( Muller et al.1 12004) g-mode oscill ations of protoneu- 



tron stars (PNSs) (|Ott et all [20061) and the sta nding- 
accretion-shock-in stability (SASI) (jKotake et all 120071 : 
iMarek et al.ll2009T ). 

However, most of them have been based on two- 
dimensional (2D) simulations that assume axisymmetry . 
Then, the growth of SASI (e.g.. iBlondin et all 12 003): 



Schcc k et all (|2004D ; lOhnishi et al.1 (|2006l ); |Fogrizzo et all 
(|2007l) ) and the large-scale convection, both of which are 
now considered to generically develop in the post bounce 
phase and to help the neutri no-driven explosions (|Bethd 
Il990t IMarek fc Janld [200l), develop along the symme- 
try axis preferentially, thus suppressing the anisotropies 
in explosions. So far very few three-dimensiona l (3D) 
studies have be en conducted (jMuller fc Jankal Il997t 
iFrver et al.ll2004 ). Moreover, neither the growth of the 
SASI nor its effects on the GWs have been studied yet 
because SASI is suppressed artificially owing to th e lim- 
ited computational do main (iMiiller fc Jankal [l997) or to 
the early shock-revival (jFrver et al. 2004) , which has not 
been discovered by other supernova simulations. 

In this Letter, we study the properties of the gravi- 
tational radiation based on the 3D simulations, which 
demo nstrate the neutrino- driven explosions aided by 
SASI (Iwakam i et al.|[20 08a) . Supported by the evidence 
of the slow rotation prior to core-collapse (jHeger et al.1 
1201)5 lQiTeTa!]l2l)0^ . we focus on the asphericities out- 
side the protoneutron stars, which are produced by the 
growth of SASI. By perfo rming a ray-tracing calculation 
in 3D (jKotake et al.l 120091 ) . we estimate accurately the 
gravitational waveforms from anisotropic neutrino emis- 
sions. We show that the features of the gravitational 
waveforms are significantly different than the ones in the 
axisymmetric cases, which should tell us the necessity of 
the 3D supernova modeling. 
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2. NUMERICAL METHODS AND MODELS 

The employed numerical methods and model con- 
cepts are essentially the same as those in our previ- 
ous paper (llwakami et all 1200 8a?). Using the ZEUS-MP 
(|Haves et all 12006ft as a hydro-solver, we solve the dy- 
namics of the standing accretion shock flows of mat- 
ter attracted by the protoneutron star and irradiated 
by neutrinos emitted from the protoneuton star. We 
employ the so-called li ght-bulb approximation (see, e.g., 
Uanka fc Mulleri |l996)) and adjust the neutrino lumi- 
nosities from the PNSs to trigger explosions. 

The computational grid is comprised of 300 logarith- 
mically spaced, radial zones to cover from the absorbing 
inner boundary of ~ 50 km to the outer boundary of 
2000 km, and 30 polar (9) and 60 azimuthal (4>) uni- 
form mesh points, which are used to c over the whole 
solid an gle (see for the resolution tests in llwakami et al.1 
(12008am . 

The initial conditi o ns ar e provided in the same man- 
ner of lOhnishi et al.l (|2006f) . which describes the spheri- 
cally symmetric steady accretion flow throu gh a stand- 
ing shock wave (jYamasaki fc Yam ada 2005). To induce 
non-spherical instability we add random velocity per- 
turbations to be less than 1 % of the unperturbed ve- 
locity. Changing the input (electron-)neutrino luminos- 
ity at the surface of the protoneutron star from L v — 
6.4 — 6.85 x 10 52 erg s _1 , we have computed six 3D mod- 
els and one 2D model (see Table [TJ) Except for model 
E, we can observe the continuous increase of the average 
shock radius with the growth of SASI, reaching the outer 
boundary of the computational domain with the explo- 
sion energy of ~ 10 51 erg. Until this moment, we run the 
simulations (e.g., At in the table). 

Following t he method in lEpsteinl (|1978f) ; 

iMuller fc Jankal (|1997l ). the two modes of the GWs 
from anisotropic neutrino emissions are derived as 
follows, 

h+ = C f f dn'(l + s(O')c((f>')s(O+c(0')c(O) x 

JO J4TT 

(fMW6 - c(fl'MQ) 2 - s 2 (flV(<//) di v (n',f ) 

[s(e')c(^)c^) - c((9')s(£)] 2 + s 2 (6')s 2 {(t>') dfl' 

and 

h x =2C f f dn'(l + s(6')c(cb')s(Z)+c(9')c(0) x 

J0 J At: 

a (g>(^)( a (fl')c(^)c(fl - c(fl') a (Q) dz„(?r , t>) 

{s(6')c(<j)i)c(£) - c(6»')s(C)] 2 + s 2 {0')s 2 (</>') dfl' ' 

where s(A) = sm(A),c(B) = cosB, C = 2G/(c 4 R) with 
G, c and R, being the gravitational constant, the speed of 
light, the distance of the source to the observer respec- 
tively, dlv/dtt represents the direction-dependent neu- 
trino luminosity emitted per unit of solid an gle into direc- 
tion o f fi, and £ is the viewing angle (e.g.. iKotake et aTl 
(2009)). For simplicity, we consider here two cases, in 
which the observer is assumed to be situated along 'po- 
lar' (£ = 0) or 'equatorial' (£ = ir/2) direction. To de- 
termine dl u /dfl, we perform a ray-tracing calculation. 
Since the regions outside the PNSs are basically thin 
to neutrinos, we solve the neutrino transport equations 
in a post-processing manner to estimate the neutrino 
anisotropy, in which neutrino absorptions and emissions 



TABLE 1 

Model Summary 



Model 




At 


\ ' 'max | 


1 ^max 1 






(10 52 erg/s) 


(ms) 


(HT 22 ) 


(io- 22 ) 


(lO- n M 


A 


6.8 


507 


3.0(x) 


3.1(X) 


6.0 


B 


6.766 


512 


2.4(+) 


3.1(x) 


7.4 


C 


6.7 


532 


3.7(x) 


2.6(x 


7.4 


D 


6.6 


667 


2.0(x) 


2-0(+) 


10.8 


E 


6.4 




0.7+) 


1.6 (x) 


2.6 


F 


6.85 


429 


2.1(x) 


4.3(+) 


5.6 


A'(2D) 


6.8 


509 




7.7(+) 


1.1 



Note. — L v denotes the input luminosity. At represents the simu- 
lation time. Model A' is the 2D version of model A. |/in lax | and l^i^ax! 
represents the maximum amplitudes seen from the pole or the equator, 
respectively. + or X indicates the polarization of the GWs at the maxi- 
mum. Eqw is the total radiated GW energy. The supernova is assumed 
to be located at a distance of 10 kpc. 

by free nucleons, dominant out s ide th e PNSs, are taken 
into account (see IKotake et al.l (|2009t ) for more details). 
To get the numerical convergence, we need to set 45000 
rays to determine the luminosity for each given direc- 
tion of f2. The GW amplitudes from mass motions 
are extrac ted by using the stan dard quadrupole for- 
mula (e.g..lMiiller fc Jankal (|1997h ) . and their spectra by 
iFlanagan fc Hughes! (| 19981 ) with the FFT techniques. 




Fig. 1. — Four snapshots of the entropy distributions of model 
A. The second and fourth quadrant of each panel shows the sur- 
face of the standing shock wave. In the first and third quadrant, 
the profiles of the high entropy bubbles (colored by red) inside the 
section cut by the ZX plane are shown. The side length of each 
plot is 1000km. The insets show the gravitational waveforms from 
anisotropic neutrino emissions, with '+' on each curves represent- 
ing the time of the snapshot. Note that the colors of the curves 
are taken to be the same as the top panel of Figure [2] 



3. RESULTS 

Figure Q] shows the 3D hydrodynamics features of SASI 
from the early phase of the non-linear regime of SASI 
(top left) until the shock break-out (bottom right) with 
the gravitational waveform from neutrinos inserted in 
each panel. 

After about 100 ms, the deformation of the standing 
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shock becomes remarkable marking the epoch when the 
S ASI enters the non- linear regime (top left of Figure Q]) . 
At the same time, the gravitational amplitudes begin 
to deviate from zero. Comparing the top two panels in 
Figure [2j which shows the total amplitudes (top) and 
the neutrino contribution only (bottom), it can be seen 
that the gross structures of the waveforms are predomi- 
nantly determined by the neutrino-originated GWs with 
the slower temporal variations (> 30 — 50 ms), to which 
the GWs from matter motions with rapid temporal vari- 
ations (< 10 ms) are superimposed. So, we first pay at- 
tention to the neutrino GWs in the following and discuss 
the importance of the matter GWs later in the spectrum 
analysis. 

As seen from the top right through bottom left to right 
panels of Figure [U the major axis of the growth of SASI 
is shown to be not aligned with the symmetric axis (:Z 
axis in the figure) and the flow inside the standing shock 
wave is not symmetric with respect to this major axis 
(see the first and third quadrant in Figure [1]). This is a 
generic feature in the computed 3D models, which is in 
contrast to the axisymmetric case. The GW amplitudes 
from SASI in 2D showed an increasing trend with time 
due to the sy mmetry axis, along which S ASI can develop 
preferentially (jKotake et al.ll2007ll2009l) . Free from such 
a restriction, a variety of the waveforms is shown to ap- 
pear (see waveforms inserted in Figure 1). Furthermore, 
the 3D standing shock can also oscillate in all directions, 
which leads to the smaller explosion anisotropy than 2D. 
With these two factors, the maximum amplitudes seen ei- 
ther from the equator or the pole becomes smaller than 
2D (compare |/imaxl m Table [T] for models A and A'). On 
the other hand, their sum in terms of the total radiated 
energy (£qw in Table Q]) are found to be almost com- 
parable between 2D and 3D models, which is likely to 
imply the energy equipartition with respect to the spa- 
tial dimensions. It was reported in a pioneering 3D study 
(|Muller fc Jankalll997f) that the emitted GW energy in 
the forms of neutrinos was reported to be about 2 or- 
ders of magnitude smaller for 3D than in 2D. However 
this may be because the global anisotropy of the SASI 
was not captured in their simulations with the limited 
computation domain in the polar direction. 

Figure [2] shows the gravitational waveforms for differ- 
ent luminosity models. The input luminosity for the two 
pair panels differs only 0.5%. Despite the slight differ- 
ence, the waveforms of each polarization are shown to ex- 
hibit no systematic similarity when seen from the pole or 
equator. This is due to the chaotic nature of SASI influ- 
enced by small differences. No systematic dependence of 
the maximum amplitudes and the radiated GW energies 
on the neutrino luminosities is found among the mod- 
els investigated here (see Table 1). In fact, the largest 
emitted GW energy is obtained for model D with the 
intermediate neutrino luminosity (Table 1). It is noted 
that the waveforms are even sensitive to the difference 
of the grid zoning due to the stochastic nature, although 
no significant difference in the GW amplitudes is seen in 
the numerical run with twice as many grid points in the 
9 and (f> directions. 

Now we move on to discuss the features of the wave- 
forms by the spectrum analysis. From Figure [3j it can 
be seen that the neutrino GWs, albeit dominant over the 
matter GWs in the lower frequencies below ~ 10 Hz, be- 
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Fig. 2. — Gravitational waveforms from neutrinos (bottom) and 
from the sum of neutrinos and matter motions (top) , seen from the 
polar axis and along the equator (indicated by 'Pole' and 'Equator') 
with polarization (+ or X modes) for models A and B. The distance 
to the SN is assumed to be 10 kpc. 
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Fig. 3. — Characteristic gravitational-wave spectra from 
anisotropic neutrino emissions ('Neutrino') and matter motions 
('Matter') with optimal detection limits of TAMA, the first LIGO, 
advanced LIGO, and LCGT for a supernova at a distance of 10 
kpc. 



come very difficult to detect for ground-based detectors 
whose sensitivity is limit ed mainly by the seismic noises 
at such lower frequenc ies (lAndo et alj2005t I Abbott et all 
120051: lWeinsteinll2002tlKuroda et al.lll999D . On the other 
hand, the gravitational-wave spectra from matter mo- 
tions peak near ~ 100 Hz, reflecting the growth of I = 2 
mode of SASI with timescales of O(10) ms. Such sig- 
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nals from a Galactic supernova are probably within the 
detection limits of the LIGO-class detectors, and seem 
surely visible to the next-generation detectors such as 
the advanced LIGO and LCGT. It is noted that another 
peak in the GW spectra near ~ 1kHz is from the rapidly 
varying matter motions of O(ms), induced by the local 
hydrodynamical instabilities. These gross properties in 
the GW spectra are found to be common to the other 
luminosity models. Thus the peak in the spectra near 
~ 100Hz is a characteristic feature obtained in the 3D 
models computed here. 

4. SUMMARY AND DISCUSSION 

We studied the properties of GWs based on the 3D 
simulations, which demonstrate the neutrino-driven ex- 
plosions with SASI. We focused on the asphericities out- 
side the protoneutron stars, which are produced by the 
growth of SASI. By performing a ray-tracing calculation 
in 3D, we estimated accurately the gravitational wave- 
forms from anisotropic neutrino emissions. Our results 
showed that the waveforms vary much more stochasti- 
cally than for 2D because the explosion anisotropies de- 
pend sensitively on the growth of the SASI which devel- 
ops chaotically in all directions. From the spectra analy- 
sis, it was pointed out that the gravitational- wave spectra 
from matter motions peak near ~ 100 Hz, reflecting the 
growth of t — 2 mode of SASI, and that such signals are 
probably visible to LIGO-class detectors for a Galactic 
source. 

Here it should be noted that the simulations high- 
lighted in this paper are only a very first step towards 
realistic 3D modeling of the supernova explosions. The 
approximations adopted in this paper, such as the re- 
placement of the PNS by the fixed inner boundary and 
the light-bulb approach with the constant neutrino lu- 
minosity, should be improved. Important ingredients for 



Ando, M., et al. 2005, Classical and Quantum Gravity, 22, 1283 
Abbott, B., et al. 2005, Phys. Rev. D, 72, 122004 
Bethe, H. A. 1990, Reviews of Modern Physics, 62, 801 
Blondin, J. M., Mezzacappa, A., & DeMarino, C. 2003, ApJ, 584, 
971 

Burrows, A. & Hayes, J. 1996, Physical Review Letters, 76, 352 
Burrows, A., Livne, E., Dessart, L., Ott, C. D., & Murphy, J. 2006, 
ApJ, 640, 878 

Cerda-Duran, P., Font, J. A., & Dimmelmeier, H. 2007, A&A, 474, 
169 

Dimmelmeier, H., Ott, C. D., Janka, H.-T., Marek, A., Miiller, E. 

2007, Physical Review Letters, 98, 251101 
Epstein, R. 1978, ApJ, 223, 1037 

Flanagan, E. E., & Hughes, S. A. 1998, Phys. Rev. D, 57, 4566 
Foglizzo, T., Galletti, P., Scheck, L., & Janka, H.-T. 2007, ApJ, 
654, 1006 

Fryer, C. L., Holz, D. E., & Hughes, S. A. 2004, ApJ, 609, 288 
Fryer, C. L. 2004, ApJ, 601, L175 

Hayes, J. O, Norman, M. L., Fiedler, R. A., Bordner, J. O., Li, 

P. S., et al. 2006, ApJS, 165, 188 
Heger, A., Woosley, S. E., & Spruit, H. C. 2005, ApJ, 626, 350 
Hirschi, R., Meynet, C, & Maeder, A. 2005, A&A, 443, 581 
Hough, J., & Rowan, S. 2005, Journal of Optics A: Pure and 

Applied Optics, 7, 257 
Iwakami, W., Kotake, K., Ohnishi, N., Yamada, S., & Sawada, K. 

2008, ApJ, 678, 1207 

Iwakami, W., Kotake, K., Ohnishi, N., Yamada, S., & Sawada, K. 

2008, arXiv:0811.0651 
Janka, H.-T., & Miiller, E. 1996, A&A, 306, 167 
Kotake, K., Yamada, S., & Sato, K. 2003, Phys. Rev. D, 68, 044023 



the GW emissions su ch as the oscillations of the PNSs 
(|Burrows et al.l 120061 ) and the enhanced neutrino emis- 
sions near the equa tors deep inside the PNS observed in 
iMarek et~a l. ( 2009), cannot be treated in principle here. 
The exploration of these phenomena will require con- 
sistent simulations covering the entire stellar core and 
starting from gravitational collapse with neutrino trans- 
fer coupled to the 3D magnetohydrodynamics. Besides, 
the eff ects of rotation on the SASI are yet to be investi- 
gated (Iwakami et al. 2008b), although the construction 
of the rotating steady accretion flows is still a major un- 
dertaking. 

Bearing these caveats in mind, the stochastic nature 
of the GWs, which is illuminated here by the 3D mod- 
els, is qualitatively new. It should be also mentioned 
that the GW spectra for the acoustic mechanism showed 
a very strong peak near ~ kHz, reflect ing the g-modes 
oscillations of the PNS (|Ott et al.ll2006f ). and are signif- 
icantly different from the waveforms discussed here for 
the SASI- aided neutrino-driven models. This means that 
the successful detection of the GW signals will supply us 
a powerful tool to probe the explosion mechanism. Re- 
cently the next-generation detectors, by which the GWs 
studied here could be surely visible for a Galactic source, 
are planned to be on line soon (2014!) (lLantzl2006lh Here 
it is an urgent task for theorists to make precise predic- 
tions of the GW signals based on more sophisticated 3D 
supernova modeling. 
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